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ABSTRACT  OF  THE  DISSERTATION 


A  Generalized  Solution  to  a  Class 
of  Printed  Circuit  Antennas 
By 

Pisti  B.  Katehi-Tseregounis 
Doctor  of  Philosophy  in  Engineering 
University  of  California,  Los  Angeles,  1984 
Professor  Nicolaos  G.  Alexopoulos,  Chair 

This  dissertation  deals  with  the  theory  and  design  of 
antennas  excited  either  by  a  microstrip  transmission  line  or 
by  a  gap  generator.  The  antennas  and  the  strip  transmission 
line  may  be  embedded  inside  or  printed  on  the  substrate.  A 
theoretical  approach  is  implemented  which  accounts  accur¬ 
ately  for  the  physical  effects  involved  including  surface 
waves.  The  Green's  function  has  been  obtained  by  synthesiz¬ 
ing  the  fields  of  Hertzian  dipoles  which  are  oriented  in 
arbitrary  directions  and  which  are  printed  on  or  embedded  in 
the  substrate.  The  method  of  solution  is  based  on  solving 
the  Pocklington  integral  equation  by  employing  the  method  of 
moments  with  proper  choice  of  expansion  and  testing  func¬ 
tions.  The  excitation  mechanism  is  taken  into  account 
effectively  by  considering  it  as  part  of  the  antenna.  The 
current  distribution  is  obtained  both  on  the  transmission 
line  and  the  printed  circuit  antennas  by  matrix  inversion. 
The  method  accounts  for  conductor  thickness  and  for 


arbitrary  substrate  parameters. 

As  an  example,  printed  strip  dipoles  excited  by  a 
transmission  line  embedded  in  the  substrate  or  by  a  voltage 
gap  generator  are  considered.  Current  distribution,  self 
impedance  and  reflection  coefficient  for  the  case  of  the 
transmission  line  excitation  as  well  as  input  impedance, 
resonant  length,  resonant  resistance  and  radiation  patterns 
for  the  case  of  the  gap  voltage  excitation  are  obtained 
for  a  variety  of  antenna  arrangements.  A  serious  amount  of 
effort  is  also  being  placed  in  evaluating  the  importance  of 
higher  order  surface  wave  modes  which  are  determined  by  the 
relative  dielectric  constant  and  the  thickness  of  the  sub¬ 
strate.  Comparison  of  the  theoretical  results  to  experimental 
data  for  the  case  of  an  electromagnetically  coupled  printed 
strip  dipole  to  a  strip  transmission  line  shows  excellent 
agreement . 


CHAPTER  1 


INTRODUCTION 

Integrated  or  printed  circuit  antennas  are  a  natural 
evolution  of  integrated  circuit  components  and  are  finding 
increased  use  in  the  microwave,  millimeter  and  far  infrared 
frequency  ranges.  Therefore,  the  development  of  antennas 
which  are  amenable  to  integration  with  other  printed  circuit 
elements  is  of  significant  technological  importance. 

The  microstrip  antenna  concept  dates  back  about  30 
years  to  work  in  the  United  States  by  Deschamps  [11  and  in 
France  by  Gutton  and  Baissinot  [2].  Shortly  thereafter, 
Lewin  [3]  investigated  radiation  from  stripline  discontinui- 
ties.  Additional  studies  were  undertaken  in  the  late  1960 's 
by  Kaloi  [54]  who  studied  basic  rectangular  and  square  con¬ 
figurations  (patches).  However,  the  inherent  advantages  of 
antenna  elements  (conformality  to  a  given  surface,  light 
weight,  negligible  volume,  inexpensiveness),  were  not  put  to 
widespread  practice  until  the  1970's  [41-121].  The  environ¬ 
mental  and  technological  constraints  having  been  resolved, 
the  task  remained  to  develop  analytical  methods  which  would 
provide  accurate  design  criteria.  Mathematical  modeling  of 
the  basic  microstrip  radiator  was  carried  out  initially 
either  by  the  application  of  transmission  line  analogies  to 
simple  rectangular  patches  fed  at  the  center  of  a  radiating 
wall,  or  by  an  open  resonator  model  [9]- [21].  The  former 
approach  gives  a  heuristic  explanation  of  the  radiation 


properties  of  the  antenna,  while  the  latter  provides  a  more 
accurate  prediction  of  the  antenna  characteristics.  However, 
both  models  apply  mainly  to  the  dominant  resonator  model  and 
their  accuracy  is  questionable  for  higher  order  modes, 
especially  because  they  do  not  account  for  the  excitation 
of  surface  waves. 

Surface  waves  have  an  important  effect  on  the  printed 
circuit  antenna  current  distribution,  as  well  as  input 
impedance,  resonant  length,  bandwidth  and  efficiency  [34]- 
(37].  In  addition,  since  surface  waves  are  cylindrical 
waves  in  nature,  they  decay  only  as  the  inverse  square  root 
of  the  distance  from  their  source  and  for  this  reason  they 
can  be  significant  in  mutual  impedance  computations  [26], 
an  Important  parameter  in  phased  array  design.  It  has  been 
shown  that,  regardless  of  the  substrate  thickness  and 
dielectric  permittivity,  the  dominant  surface  wave  mode  is 
always  excited.  The  power  propagating  in  this  mode  is  a 
function  of  the  characteristics  of  the  substrate.  As  more 
energy  is  trapped  in  the  substrate,  the  microstrip  antenna 
•ecomes  less  efficient  (35].  In  many  applications,  such  as 
in  the  millimeter  or  far  infrared  region  I29]-I31],  today's 
technology  provides  substrates  which  are  several  wavelengths 
thick.  This  permits  the  propagation  of  many  TM  and  TE  waves 
in  the  substrate,  further  complicating  the  design.  These 
■odes  can  also  cause  impairment  of  efficiency.  It  becomes 
evident  from  this  discussion  that  a  theoretical  approach 
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should  be  implemented  which  accounts  accurately  for  all  the 
physical  effects  involved  including  surface  waves.  Such  an 
approach  excludes  either  of  the  previously  mentioned  tech¬ 
niques  and  relies  on  treating  the  microstrip  element  as  an 
antenna  rather  than  a  transmission  line  section  or  as  a 
resonator. 

A  microstrip  antenna  is  usually  excited  either  by  the 
inner  conductor  of  a  coaxial  transmission  line  [39] 

(Fig.  la)  or  by  a  microstrip  transmission  line  [4l]-I43] 
printed  or  embedded  in  the  substrate  (Fig.  lb).  From  these 
two  ways  of  excitation,  the  latter  has  demonstrated  that  a 
microstrip  antenna  electromagnetically  (EM)  coupled  to  a 
microstrip  line  makes  a  desirable  element  for  one-  and  two- 
dimensional  antenna  arrays. 

Recently,  Oltman  and  Huebner  [40],  and  later.  Stern 
and  Elliott  [4l]-[42]  experimentally  studied  this  radiator 
as  an  element  as  well  as  part  of  a  two-dimensional  array  and 
they  described  a  design  procedure  with  the  objectives  of  an 
input  match  and  a  desired  radiation  pattern. 

In  this  present  work,  strip  dipoles  printed  or  embedded 
in  the  substrate  excited  either  by  a  gap  generator  or  a 
microstrip  transmission  line  are  considered.  The  thickness 
of  the  strips  is  considered  finite  and  the  widths  of  the 
dipole. and  transmission  line  are  assumed  to  be  much  smaller 
than  the  wavelength  so  that  the  transverse  components  of 
the  current  give  a  second  order  effect.  The  current  distri- 
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button  on  the  antenna  is  obtained  first  by  solving  a  two- 
dimensional  Pocklington's  Integral  Equation  [38].  The 
Green's  function  in  this  case  can  be  obtained  by  synthesiz¬ 
ing  the  fields  of  Hertzian  dipoles  which  are  oriented  in 
arbitrary  directions  and  which  are  printed  or  embedded  in 
the  substrate,  therefore  accounting  properly  for  all  the 
boundary  conditions  pertinent  to  the  problem. 

An  analytical  solution  of  the  two-dimensional 
Pocklington's  integral  equation  is  precluded  due  to  the 
immense  complexity  of  the  problem.  The  use  of  numerical 
techniques  which  discretize  the  integral  equation  and 
obtain  the  current  distribution  by  matrix  inversion  is 
necessary.  A  numerical  method  which  has  found  widespread 
and  successful  use  for  the  solution  of  Pocklington  type 
integral  equations  is  the  Method  of  Moments  [25]- [28].  For 
the  present  application,  the  Green's  function  relevant  to 
the  problem  is  given  by  Sommerfeld-type  integrals  which 
require  special  integration  techniques  when  field  and  source 
point  coincide  [38]. 

For  the  case  of  the  antenna  excited  by  a  voltage  gap 
generator,  the  input  impedance  is  defined  as  the  ratio  of 
the  applied  voltage  to  the  input  current  and  resonant  length, 
resonant  resistance,  bandwidth  and  efficiency  are  evaluated 
as  function  of  the  substrate  characteristics.  In  the  case 
of  a  dipole  EM  coupled  to  a  microstrip  line,  transmission 
line  (T.L.)  theory  is  used  to  derive  a  form  for  the  self- 


impedance.  The  application  of  T.L.  theory  becomes  possible 
by  virtue  the  fact  that  the  distance  of  the  feeding  line 
from  the  ground  has  been  kept  always  very  small  compared  to 
wavelength  in  the  dielectric  so  that  most  of  the  contribu¬ 
tion  for  the  electric  field  under  the  microstrip  line 
results  from  a  dominant  TEM-like  mode.  This  leads  to  a 
design  procedure  which,  for  a  given  substrate,  permits 
determination  of  the  length  of  the  dipole,  overlap  and  off¬ 
set  so  that  a  desired  input  match  is  achieved. 

A  serious  amount  of  effort  has  also  been  invested  in 
determining  the  effect  of  the  substrate  thickness  and 
relative  permittivity  on  the  radiation  properties  of  printed 
circuit  dipoles  (PCD's).  A  trade-off  between  substrate 
thickness  and  resonant  input  impedance,  bandwidth  and 
radiation  efficiency  is  presented  for  PTFE  glass,  quartz 
and  GaAs  substrates.  The  E-  and  H-plane  normalized  power 
patterns  are  also  examined  as  a  function  of  and  h,  and 
it  is  shown  that  even  for  thin  substrates,  multiple-beam 
radiation  can  result  for  certain  values  of  through  the 
excitation  of  surface  waves. 

Throughout  this  work,  the  cost  of  the  computer  programs 
was  kept  very  low  with  the  application  of  special  analytical 
and  numerical  techniques  for  the  evaluation  of  the  elements 
of  the.  impedance  matrix.  These  techniques  will  be  described 
since  they  are  quite  general  and  apply  to  any  kind  of 
printed  antenna. 


CHAPTER  2 

THE  GREEN’S  FUNCTION  AND  POCKLINGTON' S 
INTEGRAL  EQUATION 

2- 1 .  Derivation  of  the  Green’s  Function 

This  chapter  contains  a  development  of  the  Green’s 
Function  pertinent  to  the  problem  of  strip  or  patch  antennas 
printed  on  and/or  embedded  in  a  grounded  dielectric  sub¬ 
strate  of  thickness  h  and  relative  dielectric  constant 

In  order  to  formulate  the  Green’s  function,  two  ele¬ 
mentary  horizontal  electric  dipoles  (HED’s)  are  considered 
to  be  at  the  positions  (xj,  y’,  0)  and  (x^,  y^*  "^>5) 
shown  in  Figure  2.1.  The  assumed  time  dependence  is  e^“^ 
and  it  is  suppressed  throughout  the  dissertation.  The 
electromagnetic  field  at  any  point  due  to  these  two  dipoles 
is  the  superposition  of  the  fields  arising  from  each  one 
separately 


fi*  -  *  fi‘ 

(2.1) 

(2.2) 

with 

and  the  electromagnetic  fields 

in 

medium  (1)  and  (2)  respectively.  Maxwell’s  equations 

now 

take 

the  following  form; 

(2.3) 
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Figur#  2.1 

HED't  Printed  on  and  Embedded  in  a 
Grounded  Dielectric  Slab 
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(2.7) 

(2.8) 


where  the  superscript  i  indicates  the  medium  with  the  sub¬ 
script  V  indicating  the  source.  The  potential  function 

is  now  introduced.  will  turn  out  to  be  the  vector 
potential  function  in  medium  (i)  due  to  source  (v).  A 
substitution  of  (2.9)  into  the  curl  equation  for  gives 
the  result 

^xdj-jajXj)  «  0  .  (2.10) 

Since  ’^x(^«'^)  ■  0,  it  follows  that 

-  Ui  (2.11) 

with  the  scalar  potential  in  medium  (i)  due  to  the  source 
(v) .  The  relationship  between  and  can  be  obtained  by 
substituting  (2.9)  and  (2.11)  into  (2.3),  i.e. 


-y  5 

♦*0  V 


1  V 


^♦1 

(2.12) 


If  one  considers  the  Lorentz  condition 


^•XJ  -  jwc.y^jt^ 


(2.13) 
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equation  (2.12)  reduces  to 

0?^  *  (2.14) 

vdiere 

X- 

The  Lorentz  condition  enables  one  to  write  the  electromag¬ 
netic  field  in  medium  (i)  due  to  source  (v) ,  since  it 
follows  that 


l^eiVo 


IV  V 


(2.15) 


where 


<D^e  V 
2  0  0 

r  o  o 


i  -  2 


(2.16) 


For  the  case  of  infinitesimally  small  dipole  sources 
oriented  along  the  x  or  y  axis  the  solution  to  the  inhomo¬ 
geneous  equation  (2.14),  when  that  solution  also  satisfies 
the  specified  boundary  conditions,  is  called  Green's 
function  and  is  given  by 

^1 


For  an  elementary  dipole  of  arbitrary  orientation,  this 
vector  function  becomes  dyadic  and  is  the  solution  of  the 
following  inhomogeneous  equation 

^  t  j  wu 

ir  *  kpTi  -  — I  6(r  -  ?  )  (2.18) 

1  V  kj 

where  T  is  the  unit  dyadic  or  idemfactor  given  by  xx  ♦  yy+  k. 
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Considering  all  the  equations  above,  one  can  prove  145]  that 
for  the  case  of  a  general  current  the  electric  field  is 
given  by 

eJcO  -  J  th?  T  ♦  Wl  •  Pjcr/r’)  •  3^(r')dT’ 

(2.19) 

with  the  integration  extended  over  the  volume  which  contains 
all  the  current  sources.  One  can  also  deal  directly  with 
the  equation 

^2.20) 

which  the  electric  field  satisfies.  In  this  case,  the 
the  dyadic  Green’s  function  can  be  defined  as  that  solution 
to  the  inhomogeneous  equation 

-  -ju>uj6Cr-r^)  (2.21) 

which  also  satisfies  the  specified  boundary  conditions. 

Using  equations  (2.20)  and  (2.21),  one  can  prove  145],  [46] 
that 

tj,(r)  -  Gj(r/r')  •  :5^(r’)dT’  (2.22) 


Therefore,  it  can  be  concluded  that 


^i(r/r’)  -  (k?T  ♦  •  Fj(r/r') 

V?,r*eV 


(2.24) 


For  the  total  field,  equations  (2.19)  and  (2.22)  can  be 


written  as 


ti(r)  -  L  /  Cjcr/r')  •  5  (r')dT*  - 
^  v-1,2  ^  ^ 

'  E  /  (k-T  +  •  Pjcr/r')  • 

..^1  O  •'tr  *  V 


v-1,2  -V, 


(r')dT* 


(2.25) 


The  solution  which  satisfies  equations  (2.18)  and  the 
appropriate  boundary  conditions  in  medium  (i)  consists  of 
two  parts;  the  secondary  solution  which  is  the  solution  to 
the  corresponding  homogeneous  equation  and  the  primary 
solution  which  is  the  particular  solution  to  the  wave 
equation.  For  the  case  of  a  horizontal  dipole  current  along 
the  X  direction  3  (r’),  the  dyadic  function  F^(r/r’)  has  the 


form 


fiC/r')  .  XX  .  XX 


(2.26) 


where  .the  components  F^,._„(r/r’) ,  F^,-„(r/r')  were  found  to  be 


C33J.  1371 
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(2.27) 


•  {1  -  fii2*v2  *  *  u^sinhCuz'’)  1) 


and 


sinhluCh-z'^  + 


XdX 


fjCX.Ej.,  h) 

•  /juy  \  r* 


-V  *n 


>4»k^ 

sinh[u(h‘“z'^  ♦ 
f  (^ » 

cosh[u(h*z^  ♦  2^6.,)]  7 

- r-T! - ^rr-^ —  ^ 

f2{X,e  ,h) 


(2.28) 


The  various  parameters  involved  in  equations  (2.27)  and 
(2.28)  are  defined  as 


u 

*  ,  u^  »  [X^-kJ)^/^ 

(2.29) 

p 

-  [(x-x’)2  ♦  (y-y’)^]^^^ 

(2.30) 

^^(^.eT.h) 

*  u^sinh(uh)  ♦  ucosh(uh) 

(2.31) 

^2 

■  E^u^cosh(uh)  ♦  usinh(uh) 

(2.32) 

The  zeros  of  fj(X,ej.,h)  and  f2(X,Ey,h)  lead  to  TE  and 
TM  surface  wave  inodes  respectively  135),  137)-I391. 


2-2  Pocklington* s  Integral  Equation 

In  this  work,  the  thickness  of  the  metallic  strips  is 
considered  finite  and  the  widths  of  the  dipole  and  the 
transmission  line.  Fig. (2.2)  and  (2.3)  are  assumed  to  be 
■uch  smaller  than  the  wavelength  in  the  dielectric  so  that 
the  transverse  components  of  the  current  are  a  second  order 
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Figure  2.3 

Strip  Dipole  Excited  by  a  Qap 
Generator 


effect.  If  it  is  assumed  that  the  current  density  of  the 
dipole  is  given  by 


(2.33) 


where  z'’  ■  0  when  v  ■  1  and  z'^*  -bg  when  v  »  2,  then  the 
integral  equation  takes  the  form 

^^Cr)  -  L  f  J  <5y’  f  'dxWk?  T  +  • 

v-1,2  J-v^/2  Jo 


•  Fj(r/r')*xJ^(r’) 


(2.34) 


If  equation  (2.26)  is  substituted  into  (2.34)  the  expression 
for  the  electric  field  takes  the  form 


f  v^/2  rv 

r  /  V,  J 

-1,2  J-\i^/2  Jq 


dx’  . 


iH^v, 


/8^F^ 

J  vxx  . 

^Ozsr  * 


a^F^. 


J^(x’,y’,z'’)2 


(2.35) 


From  (2.35)  it  is  observed  tha^  the  electric  field  consists 

of  three  components,  namely,  E^,  E^,  and  E^.  However,  only 

X  y  z 

the  E  component  is  needed  for  the  application  of  the  method 
of  Moments  along  the  x-axis.  From  equation  (2. 35) ,  the  E^ 


component  is  given  by 
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E^C?) 


r  m^/2  fK 

^  /  V,  j 

vl,2  •iwV2  •'o 

l« 2pi  vzx  .  vz 

[^i^vxx  “^3^ 


F  I 

vzx 

Xdz 


J^(x\y’,z'') 


(2.36) 


A  consideration  of  the  following  relationship 


(2.37) 


leads  to 


Ei(?) 


^  rJri^/2  f\ 

L  /  V,  <*y’  /  dx’ 

v«l,2  J-w'’/2  ■'o 


^i^ixx  *  ^  ^^ixx'^iz^  J^(x\y’,z'') 


with  F^^  given  by  (see  Appendix  A) 


i  /3“V  \  r* 

C^ilUoCOsh(uh)-6j^2'*^^"^l^(^"^^^  0 


■  2  f  j 


*6^jSinh(uh)  ♦ 


where 


1  i  -  1 

0  i  »<  1 


1  V-  1 
0  1 


(Z.38) 


XdX  (2.39) 


(2.40) 


CHAPTER  3 


TRANSFORMATION  OF  THE  INTEGRAL  EQUATION 
INTO  A  MATRIX  EQUATION 

L .  Method  of  Moments 

The  purpose  of  this  chapter  is  to  present  the  basic 
nathematical  techniques  for  reducing  functional  equations  to 
matrix  equations.  These  techniques  are  then  applied  to  the 
specific  problem  of  the  strip  dipole  printed  on  or  embedded 
in  a  grounded  dielectric  substrate  and  excited  either  by  a 
microstrip  transmission  line  or  by  a  gap  generator.  A 
unifying  principle  for  such  techniques  is  found  in  the 
general  method  of  moments,  in  terms  of  which  most  specific 
solutions  can  be  interpreted. 

Throughout  this  chapter,  the  width  of  the  strip  dipole 
and  transmission  line  is  assumed  to  be  small  enough  compared 
to  the  wavelength  in  the  dielectric,  so  that  the  transverse 
component  of  the  current  may  be  assumed  to  be  a  second  order 
effect.  It  is  to  be  emphasized  that  in  this  dissertation 
the  thickness  of  the  conducting  strips  is  finite.  The 
method  of  moments  is  a  general  procedure  to  solve  linear 
inhomogeneous  equations  of  the  type  [47] -[52] 

Lop(f)  •  g  (3.1) 

where  is  a  linear  operator,  g  is  the  source  or  excita¬ 
tion  (known  function)  and  f  is  the  current  or  response 
(unknown  function  to  be  determined) .  The  term  deterministic 
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means  that  the  solution  to  equation  (3.1)  is  unique;  that 
is,  only  one  f  is  associated  with  a  given  g.  The  integral 
equation  (2.37)  for  the  specific  problem  that  is  studied 
here  can  be  transferred  into  an  operator  equation  as 


follows : 


L^(3)  -  eJx 


(3.2) 


with  given  by 


v«l,2  •*w'^/2  •'o  I  vj 


vxx  vz^i 


(3.3) 


In  addition  to  the  above,  it  is  necessary  to  define  the 
inner  product  <5, (E^  x)>,  which  is  a  scalar,  to  satisfy 
the  following  relations  in  Hilbert  space 


<:!,  (eJ  X-)  -  <(eJ  i),J> 

<a^  ♦  6(E^  x),h>  ■  a<J,]K>  ♦  CE^  x)  ,^> 

<?*  d>  >0  if  J  ^  0 
»0  if  5  ■  0 


(3.4) 

(3.5) 


(3.6) 


where  a, 6  are  scalars,  Jxli  ■  0  and  *  denotes  a  complex  con¬ 


jugate. 


A  suitable  inner  product  for  this  problem  is 

<5,(eJJ  x)>  -  J  3  •  (EjJ  x)dx 


(3.7) 


Furthermore,  G^(r/r')  -  G^(r'/r)  (v«l,2)  and  using  equation 
(2.25)  it  can  be  shown  that 


(3.8) 


<L^(J),  (E^  x)>  «  <J,  L^CEj  X) 

which  means  that  the  integral  operator  is  self-adjoint. 
Consideration  of  equation  (3.8)  and  of  the  fact  that  zero 
excitation  gives  no  response,  it  can  be  proved  that  there 
exists  a  unique  solution  to  functional  equation  (3.2)  and, 

•  A 

therefore, the  inverse  operator  (L^)”  exists  such  that 

J  -  (L^)’^  (Ej  X)  (3.9) 

In  order  to  obtain  a  solution  for  equation  (3.2)  in  the 
form  (3.9),  we  have  to  follow  the  procedure  described 
below: 

1)  Expand  the  unknown  vector  5  in  a  series 
of  basis  functions  spanning  J  in  the 
domain  of  L^. 

2)  Determine  a  suitable  inner  product  and 
define  a  set  of  weighting  functions. 

3)  Consider  the  inner  products  of  these 
functions  with  both  sides  of  functional 
equation  (3.2)  and  transform  it  into  a 
matrix  equation. 

4)  Solve  the  matrix  equation  for  the  unknown 
vector  5. 

3-2.  Galerkin’s  Method 

For  the  evaluation  of  the  current  distribution  on 
printed  circuit  antennas,  a  specialization  of  the  general 
method  of  -moments  is  particularly  convenient.  At  first,  the 
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unknown  current  is  expanded  in  a  series  of  functions  in  the 
donain  as  follows 


3(7)  -  £  3^?) 


v*l ,  2 
ncN 


(3.10) 


where  the  I^'s  are  constants.  The  functions  are 

n  n 

called  expansion  or  basis  functions.  For  exact  solutions, 
(3.10)  is  an  infinite  summation  while  for  approximate  solu¬ 
tions  it  is  usually  a  finite  summation.  Substituting  (3.10) 
into  (3.2)  and  using  the  linearity  of  L^,  one- can  have 

L  L^^nCr))  -  (Ej  X)  (3.11) 

n 

Furthermore,  a  set  of  weighting  functions,  or  testing 
functions,  is  chosen  to  be  identical  with  the  basis 
functions,  i.e.. 


{wj^(r)}  »  {5^(r)}  m,  neN 

and  then  the  inner  product 


(3.12) 


(3.13) 


is  formulated.  This  set  of  equations  can  be  written  in 
matrix  form  as 


where 


('nl  - 

(3.14) 

'-*(3il)> 

mn  m  n 

(3.15) 

tv^i  -  <«i,  (e‘  J)> 

(3.16) 
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Since  the  matrix  (as  it  will  be  shown  in  the  follow- 

mn 

ing  chapters)  is  nonsingular,  its  inverse  exists  and 
equation  (3.14)  gives 

tin'  ■  •''b'  (5-II' 

The  vector  is  written  as 

-  tJj.  Jj*--**  ^1»  (3.18) 

and,  therefore,  the  solution  to  equation  (3.2)  is  given  by 

3(?)  - 13;;)'^  [z^r^  [vii  ts.is) 

One  of  the  main  tasks  in  any  specific  problem  is  the  choice 

of  the  functions  5^.  They  should  be  linearly  independent 

and  selected  so  that  some  combination  (3.10)  can  approximate 

J(r)  reasonably  well.  Additional  factors  which  affect  the 

choice  of  are: 
n 

i)  The  desired  accuracy  of  the  solution, 

ii)  The  ease  of  evaluation  of  the  matrix  elements, 

iii)  The  size  of  the  matrix  that  can  be  inverted, 

iv)  The  convergence  of  the  solution. 

3-3  Impedance  Matrix  Element  Formulation 

As  mentioned  previously,  the  choice  of  the  basis 
functions  is  determined  by  many  factors  dictated  by  the  pro¬ 
blem  under  consideration.  For  the  case  of  a  printed  or 
mmbedded  strip  dipole  electromagnetically  coupled  to  a 
■icrostrip  transmission  line  or  excited  by  a  gap  generator 
■the  basis  functions  were  chosen  to  he  of  the  form  _ _ 

22 


(3.20) 


“  xjJJCr*)  -  xJ^3j(x’)J^(y)6(2’-2'’) 

The  overlapping  piecewise  sinusoidal  functions 

(381  of  the  form  (Fig.  3.1), 


sinIk(x’-Xj^_j)]  ^  sin[k(Xj^^j-x')l 
^n-1  sin(M  )  *  ^n  sin(let^) 

^  X 


(3.21) 


with 


1  X„  1  <  x*  <  x^ 
n-1  —  —  n 

0  elsewhere 


1  «„  i  1 

0  elsewhere 


(3.22) 


aeR'*’ 


(3.23) 


Jy(y’)  gives  the  correct  transverse  variation  of  the 
current  density  on  the  strips  taking  into  account  finite 
conductor  thickness.  The  expression  for  the  function 
Jy(y')  is  given  by  (Fig.  3.2) 


J^Cy’) 


^  -(f/l 


(3.24) 


Here,  is  the  effective  strip  width  given  by  *  w'’  +  26^’. 
The  parameter  5^  is  the  excess  half  width,  and  it  accounts 
for  fringing  effects  due  to  conductor  thickness  (53]. 
Interpretation  of  the  choice  ^or  the  current  density 
dependence  in  y* ,  indicates  that  the  edge  condition  is 
satisfied  at  y*  ■  ±W^/2,  which  is  an  equivalent  strip  of 
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Figure  3.1 


Piecewiee-Sinueoidal  Currents  on  the  Printed 
Dipole  and  the  Embedded  Microstrip  Line 


I 


/  /  /  /V  /  /  y  /  /  / 


Figur#  3.2 


Current  Dittributfon  on  tho  Dipolo  in  the 
Tranovoroo  Direction 


zero  thickness.  At  y*  ■  the  current  density  remains 

finite,  as  is  the  case  for  nonzero  thickness  conductors. 

This  choice  of  dependence  in  y*  for  the  current  density 
yields,  as  will  be  shown,  very  accurate  results  for  micro- 
strip  dipole  resonant  length. 

If  one  substitutes  (3.20)  and  (3.10)  into  (2.38),  then 
it  can  be  found  that  the  expression  for  the  electric  field 
is 


♦  Fj^^(x-x’,x^)l  ♦  k2lFjjx+x',x^)  ♦ 


♦  fJj(x-x’ ,Xj^)l  sin[le(tjj-x')l  + 


*  ISCx’+tjj)  + 


+  fi(x'-£jj)  -  2cos(k£^)6(x’)l 


(3.25) 


Now,  the  following  inner  produce  is  defined 


r*m 

x)>  -  J  J^(x)  Ej(r)6(y-y^)6(z-z^)  (3. 


Equation  (3.26)  combined  with  (3.25)  and  (3.15)  gives  the 
following  form  for  the  elements  of  the  impedance  matrix 


’»’  tsii>(ktj,)l^  /2  1  j 

*  © 

/I 

*  dx  j  dx'  (kf-k2)  IFi^^(x«x',x„.x„)  * 
u  O 


*  ''vxxf*-*'*S-V  *  ''vxx(-='**'-'‘n-V  * 

*  FijCx-x'.Xn.V  "  •’'n-V’  * 

+  sinik (ijj-x)  1  •  sinllf (ijj-x*)  1 

TL  ^  •“  I 

r^x 

’'J  *  Kxx^'^**'  ’ 

o 

-  4 ’*nV  -  '=«'■***■  •■‘n-’'.;”  ■ 


with 


•  sin(k(t^-x)l  t6(x’  +  0  ♦  6(x’-l  )  - 


2cos  (kfcj^)6  (x’ )  1 


(3.27) 


“  ^vxx  lp.{[x+x’)+(x^-Xj^)]^+(y-y’)^} 


(3.28) 


■  I  7  2  ^ 

'p  -  {(x*x*)  +  (x„-xjr+(y-y’)  ) 


(3.29) 
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Figure  3.3 

Printed  Strip  Dipoie  EM  Coupied  to  an  Embedded 
Microetrip  Line  which  ie  Fed  by  a  Coaxial  T.L. 


A  mentioned  in  previous  chapters,  the  strip  dipole  is 
excited  either  by  a  voltage  gap  generator  or  by  a  nicrostrip 
transmission  line  embedded  in  the  dielectric.  In  the  first 
case,  the  gap  generator  is  considered  at  the  feed  point  and 
all  the  elements  of  the  excitation  vector  become  zero 
except  the  one  which  corresponds  to  the  infinitesimally 
small  gap.  In  the  case  of  the  electromagnetically  coupled 
strip  dipole,  the  problem  of  the  excitation  of  the  trans¬ 
mission  lint  has  to  be  resolved.  In  practice,  the  micro¬ 
strip  line  is  excited  by  the  inner  conductor  of  a  grounded 
coaxial  line  as  shown  in  figure  (3.3).  Since  the 
theoretical  analysis  for  such  an  excitation  is  very 
difficult,  other  possible  models  were  studied  such  as  a 
voltage  or  current  generator  at  the  end  of  the  microstrip 
or  a  voltage  gap  generator  two  subsections  from  this  end. 
From  these  models  the  only  one  that  gave  very  good  results 
was  the  gap  generator  but  with  the  condition  that  the  length 
of  the  microstrip  line  is  chosen  to  be  more  than  three  wave¬ 
lengths  in  the  dielectric.  This  model  has  been  used  for  the 
derivation  of  the  results  which  will  be  presented  in  follow¬ 
ing  chapters. 


CHAPTER  4 


EVALUATION  OF  THE  SOMMERFELD 
TYPE  INTEGRALS 

Singular  Points  and  Related  Surface  Waves 

As  shown  in  Chapter  3,  the  elements  of  the  generalized 


impedance  matrix  are  given  by 


where^^^  is  a  multiple  space  integral  operator  acting  on 
the  zeroth  order  Bessel  function,  A^'^(X,ej.,h;z^,z'^)  is  a 
complicated  expression  of  transcendental  functions  without 
singularities  and  fj^(X,Ej.,h),  f2(^»Ej.,h)  are  given  by  (2.30) 
and  (2.31).  The  integral  in  equation  (4.1)  is  a  Sommerfeld 
type  [44]  and  the  existence  of  essential  singularities  in 
its  integrand  necessitates  very  careful  treatment.  In  this 
chapter,  the  computation  of  this  integral  will  be  shown 
explicitly  and  the  approximations  employed  will  be  justified 
with  estimation  of  the  introduced  error. 


In  equation  (4.1),  the  integrand  is  a  function  of 
parameter  X  through  the  radicals 


u  -  [X^-k^]’* 

(4.2) 

u©  -  tX^-kJ]** 

(4.3) 
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where  1^2  are  the  vacuum  and  substrate  wavenumbers 
respectively.  The  sign  of  the  radical  u  does  not  affect  the 
single  value  of  the  integrals,  as  the  terms  involving  the 
radicals  are  even  functions  of  u.  For  this  reason,  only 
the  branch  cut  contribution  by  the  radical  u^^  is  considered 
and  its  direction  is  determined  by  the  requirement  that  the 
radiating  field  is  a  wave  receding  from  the  source.  As  a 


result,  the  restrictions  on  X  are: 

Re(X)  >  0  (4.4) 

Im(X)  >  0  (4.5) 

which  in  turn  impose  the  following  behavior  for  u  and  u^: 

Re(u^)  >  0  Iro(u^)  >  0  (4.6) 

Re(u)  >  0  Im(u)  >0  (4.7) 


A  possible  position  of  the  branch  cuts  governed  by  these 
inequalities  is  shown  in  Fig.  (4.1). 

The  integrand  in  (4.1)  has  poles  whenever  either  one 
of  the  functions  f^(X,Cj.,h),  f2(X,Ej.,h)  becomes  zero.  The 
zeros  of  these  two  functions  correspond  to  surface-wave 
nodes.  Particularly  the  zeros  of  fj^(X,Ej.,h)  correspond  to 
TE  surface  waves  and  the  zeros  of  f2(^»Ey,h)  to  TM  surface 
waves.  In  the  case  of  a  lossless  dielectric,  these  TE  and 
TM  poles  are  the  roots  of  the  equations  u^  ■  ucoth(uh)  and 
■  •utanh(uh).  Furthermore,  these  poles  lie  within  the 
range  kj  <  Re(X)  ^  k2. 
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The  semi-infinite  integration  in  equation  (4.1)  is 
performed  along  the  real  axis  and  is  completed  in  two  steps. 

(i)  Numerical  integration  over  the  interval  tO»A] 
where  A  satisfies  the  relationship  cothl  (A^-k^)**! 

-  1. 

(ii)  Combination  of  numerical  and  analytical  integra¬ 
tion  for  the  evaluation  of  the  tail  contribution 
which  is  actually  the  integration  over  the  path 
ra,«). 

Subsequently,  it  is  concluded  that  the  (m,n)/(i,v)  element 
of  the  impedance  matrix  can  be  split  into  two  parts,  viz.. 


Z 


iv 

mn 


mn 


(4.8) 


where 


and 


r  £  itV-op)) 

Jo  v-1,2  mn  ° 


r  r  x''  {j„(Ap)} 

Ja  v*1,2  mn 


A^''(X,ey,h;z\z'’) 


dX 


(4.9) 


A^''(X,Er,h;z\z'') 


f .  (X,ej.,h)  •f2(X,ej.,h) 


dX 


(4.10) 


4-2.  Numerical  Integration  Xe[0.A] 

The  first  part  of  each  element  of  the  generalized 
impedance  matrix  is  given  by 

L  {J^CXp))  dX 
mn 

(4.113 


/■A  ^  A^'’(X.er.h;z^,z'') 

-  X  v-^Z  [fiU,c^;h).f2(X,c^7T.7 
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and  is  evaluated  numerically.  The  operator^*  is  of  the 

mn 

form 


mn 


L 


dy’JJIcy) 


72 


L 


If- 


•  Id(p-Pi)  ♦  fi(p-p2)  ♦  4(P‘P3)  ♦  fiCp-P^)]  ♦ 

•  +  6(x’-tp  -  2cos(kt^)6(x’)l  . 

•  IfiCp-Pj)  ♦  dCp-Pj)!  (4.12) 

where 


Pi  ‘ 

{(x+x‘+x^-x^)^+  (y-y)^> 

(4.13) 

•o 

1 

1 

{(x-x’+Xjjj-Xj^)^*  (y-y’)^} 

(4.14) 

P3  ‘ 

{(-x*x’+Xjjj-Xj^)^+(y-y’)^>** 

(4.15) 

P4  • 

{ (-x-x*+Xjjj-Xj^)^*(y-y’)  > 

(4.16) 

belong  to  the  set  of  basis  functions 
which  have  been  used  for  the  application  of  the  moments 
method.  In  the  integral  of  equation  (4.11),  it  has  been 
found  that  the  space  integral  operator  can  make  the  inte¬ 
grand  a  slowly  varying  function  of  X  thereby  minimizing 
the  error  of  integration.  Therefore,  for  the  evaluation  of 
^mn*  function 

S^(X)  •£!  {Jj,(Xp)}  (4.17) 

mn 

is  transformed  into  a  rapidly  converging  series  of  the  form 
(Appendix  B) 


77  J  CXX  )  + 

o  o  mn 


■"  ^o^^^in(n-l)^ 


(4.18) 


where  i„  is  the  subsection  length  (t„  <  0.02X^),  J  (Xx  )  is 
X  X  o  o  mn 

an  expression  involving  Bessel  functions  and  is  given  by 


X  „  *  (x  -x„) 
mn  '■  m  n 


(4.19) 


For  the  integration  with  respect  to  X,  because  of  the 
existence  of  the  poles  in  the  strip  Re(X)  £  k2, 

a  further  division  of  the  integration  integral  into  sub¬ 
intervals  takes  place  as  is  shown  below: 

i)  0  £  X  £  kj,  the  integration  with  respect  to  X 
is  performed  numerically  using  a  Gaussian- 
Quadrature  method  with  fixed  ints. 

ii)  1  X  <_  k2;  for  the  integration  along  this  inter¬ 
val  a  singularity  extraction  technique  is  used 
[37] -[381  which  transforms  the  integral  into  a 
finite  series  plus  an  integral  of  a  slowly  vary¬ 
ing  function.  This  finite  series  gives  the  con¬ 
tribution  of  the  surface  wave  modes  with  the  num¬ 
ber  of  its  terms  dependent  on  the  thickness  of  the 
dielectric  as  well  as  the  dielectric  constant, 

iii)  1^2  —  ^  integration  is  invoked  here 
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in  exactly  the  same  manner  as  it  is  performed  in  the  first 
sub interval . 


4-3.  Evaluation  of  the  Tail  Contribution 


In  this  case,  the  integration  with  respect  to  X  is 
extended  along  the  interval  lA,*#] .  The  integrals  which  have 
to  be  computed  are  given  by 

i  r*  r-  A^'’(X.e  ,h;2^,z'')  v 

■  i  v.f.2 


(4.20) 


where 


(kJ-k2)+E2l?]  • 

•  ^l-^i2*v2  '*■  *  UpSinhCuz'’)]^* 

•  sinh[u(h-z'’  +  * 

+  [Ejlf^-E2k]  •  |6j2U^cosh(uh)-6^2'*®^’''^I'^(^"^^)ij  ’ 

•  (Cj.-1)  •  j^6^jSinh(uh)  +  6^2sinh[u (h-z'^)  ij  j  X 


(4.21) 


with  Ej,  E^  constants  which  take  the  values  1,0.  When  the 
source  point  and  the  observation  point  coincide,  p-^0,  z^  » 
and  the  integral  above  takes  the  form. 


zi^c-)  •  c 


«  A^''(X,ey.h;z^.z'') 

Ja  v»1,2 

m  « 


(4.22) 
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where  C  is  a  constant.  For  the  integrand  in  (4.22),  it  can 


be  proved  that,  when  A  <<  X, 

A^'’(X,e^,h;2\2'') 
fjt^.eyfhj •t2lX,e^hJ  “ 


ElQlCX.ey 


,h)e 


+E2Q2(^.erb)e 


£j(X,Cj.,h)  •  f2(X,ej.,h) 


Ej^Qj^  (X  ,c^,h) *E2Q2 CX , E^,h) 

fj^CXyC^jh)  •f2(X,ej,,h) 


0(X"®)  (a<l) 


(4.23) 


where  Qj^(X,Cy,h)  and  Q2(^»Ey»b)  are  nonsingular  functions 
of  X,Cj.,h  as  a  consequence  of  which  becomes  infinite.  In 
order  to  avoid  this  difficulty,  which  arises  in  the 

computation  of  the  diagonal  elements  of  the  impedance 
matrix,  the  conducting  strips  are  assumed  to  have  a  very 
small  but  finite  thiclcness,  t,  with  the  source  current 
flowing  on  the  bottom  surface  and  the  observation  points 
located  on  the  top  surface.  As  a  result,  the  distance 
between  source  and  observation  points  is  prevented  from 
becoming  tero;  the  minimum  value  it  can  take  is  equal  to 
the  thickness  of  the  strip.  Since  A  has  already  been 
chosen  in  such  a  way  that 

cothI(A^-k2)‘‘h)  •  1  (4.24) 
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equation  (4. 20) can  be  approximated  as  (Appendix  C) , 


EiV( 

am' 


« 

E  D,(A)  i 

V4irkr/ 

v-1.2  1  ^  J 

mn 


•  du„  ♦  D,(A) 


mn 


-uC2b.-t)  V 

e  ^  {JjjCupp}  du  + 


mn 


♦  Error 


(4.25) 


where 

Dj(A) 


llEjCkJ-V^)  .  EjVl 

I 

•  e;;:^  •  l-e^lA))! 

BzW  •  j  [EjCkJ-k^J.EjV]  ^ 

+  lEjk^-Ejk]  (*ii*v2'‘*i2®vl^(r*e^TO'  '  ' 

•  rVe^(A))[  (‘-27) 

|-iEj(k|-k2)  .  Ejki  *  «i 


B.CA) 


'i2®v2 


and 


(iTe^  -  ■  rreTwr)! 

.2 


e,(A)  - 


1  *^1 


(4.28) 


(4.29) 
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(4.30) 


6  2  (A) 


ejCA) 


64  (A) 


1  *^2  *^1 


A^-k! 


k^-k^ 
*2  *1 


A^-k 


k^-k^ 
1  *^2  *1 


A^-kt 


(4.31) 


(4.32) 


er(A) 


CfiCA) 


A^-kt 


k^-k^ 

K2  Kj 


(4.33) 


(4.34) 


pj  •  ptl+ej(A)l 
pj  ■  pri+e5(A)] 
t*  -  bg  +  tll+2e4(A)) 


(4.35) 


(4.36) 


(4.37) 


The  error  made  by  the  approximation  in  equation  (4.25) 


depends  on  A.  It  has  been  found  that 


Error  -  0(A’2|z"1) 


(4.38) 


and  for  this  reason  A  always  takes  values  larger  than  10^. 
Therefore,  the  approximation  considered  here  is  a  very  good 
one  since  the  overall  error  made  in  the  computation  of 
the  input  impedance  is  of  the  order  of  0.11.  In  equation 
(4.25),  the  integrals  can  be  put  in  the  general  form, 

I.-  f*2  2  fc  ®  {Jo(«iPn}<iu.  (i-1,2)  (4.39) 

^  J(A^-kr)’»  mn  °  ^  ^  1 


where  ■  (X^-kp^  ;  a  can  take  either  one  of  the  posi¬ 
tive  values  t,  t*,  2bg-t.  If  one  interchanges  the  order  of 
the  operator  and  the  integration,  the  integral  can  be 
written  as 


'i 


£ 

mn 


mn 


1701  * 

(^o/g^ 


-r(:— 

•*0  \  gj 


where 


g^fAl 


1  ♦ 


^i 


(A) 

h 


0  V 

P  + 

gTOTj 

(A))> 

> 

I'-t 

mn 

>  - 


(4.40) 


(i  '1,2) 


(4.41) 


If  one  substitutes  (4.40)  into  (4.25),  2^(*)  becomes 


Z^'’(-) 

mn 


£ 

v-1,2 


D^A)  V 

il^^mn 


!(•■•  NT 


D^tA)  V 

D,(A)  ,v 


*2^  nn 


^  fA 

*  £  /  Dj(A)  \ - g-nn — 

vl,2  Jo  ^  ^ 


|(•'  -[i^] ) 

IP 


I  * 


w 


-) 
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♦  DjCA) 


-X(2b5-t)/g2(A) 


g2CA) 


) 


L  {J^ap)}  d\ 

inn 


(4.42) 


From  equations  (4,8),  (4.11)  and  (4.42),  it  can  be  seen 
that  the  elements  of  the  impedance  matrix  can  be  split 
into  two  other  parts  so  that  the  evaluation  of  each  is 
simplified 


2“  -  ^  (A)  (-) 


(4.43) 


mn 


mn 


where 


(A) 


mn 


r 


I 

v-1,2 


B^''(X,e^,h;2\z'’) 


L  {J  (Xp)>dX 
mn 

(4.44) 


^  (*)  ■ 

L  , 

' 

Di(A)  ,v| 

'■f 

mn 

v-1,2  " 

l'’  [g^wj 

) 

D2(A)  ,  V 

mn 


D3(A)  ,  V  2 


V 

mn 


l&S]’) 


(4.45) 


and 


B^'’(X,Cy,h;z^,z'') 


.iv 


A*''(X.Cy,h;z^z'') 
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The  first  term  on  the  RHS  of  equation  (4.43)  is  evaluated 
with  the  technique  which  was  described  in  Part  4.2  of  this 
chapter,  while  for  the  evaluation  of  the  second  term  a  more 
careful  treatment  is  necessary  in  order  to  keep  the  error 
of  integration  within  acceptable  limits.  The  operation 
of^  on  the  function 
written  in  the  final  form 


where  a^  -  constant,  ♦  *m’*n  g^T^r) 

a  slowly  varying  function  of  p.  In  equation  (4.47)  all  the 
remaining  integrals  are  evaluated  numerically  by  using  the 
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I 


Gaussian  Quadrature  method  of  integration  and  in  this 
manner  the  error  introduced  is  of  the  order  of 

'  mn  / 


I 


CHAPTER  5 


NUMERICAL  RESULTS  FOR  SELF  AND  INPUT 
IMPEDANCE  OF  PRINTED  STRIP  DIPOLES 

5-1.  A  Printed  Strip  Dipole  Excited  by  a  Microstrip 

Transmission  Line  Embedded  in  the  bielectric 

In  this  case,  the  excitation  mechanism  is  provided  by 
a  strip  transmission  line  embedded  inside  the  substrate  (see 
Figures  5.1,  5.2)  which  couples  energy  parasitically  to  the 
microstrip  dipole.  Since  the  radiation  mechanism  of  a 
microstrip  dipole  is  very  similar  to  that  of  a  microstrip 
patch  [55],  the  model  developed  in  this  dissertation  is 
applicable  to  the  analysis  and  design  of  microstrip  elements 
which  are  rectangularly  shaped  but  with  a  width  smaller  than 
the  element  length.  Reference  to  Figures  5.1  and  5.2  gives 
reinforcement  to  the  assertion  that  the  parameters  of  the 
problem  are  arbitrary  in  the  development  of  the  model, 
including  substrate  thickness  and  relative  permittivity  to 
account  for  previously  recognized  substrate  effects  [33]  - 
[38].  In  addition,  the  thickness  of  the  metallic  conductors 
is  included,  strip  transmission  line  and  microstrip  antenna 
widths  may  differ  and  the  effects  of  the  microstrip  dipole 
overlap  and  offset  with  respect  to  the  transmission  line  on 
the  current  distribution  are  investigated.  If  one  deter¬ 
mines  the  current  distribution  by  applying  the  method  of 
moments  (see  Chapter  3),  then  transmission  line  theory  is 
invoked  to  evaluate  at  the  chosen  reference  plane,  the  self 
impedance  of  the  microstrip  dipole.  This  leads  to  a 
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Figure  6.1 


8idg  and  Top  Viow  of  o  Printod  Strip  Dipolo 
Exeltod  by  a  Tranomloslon  Lino  Emboddod  In 
-  ttio  DIoloetrie  _ 


offset 


Figure  5.2 

Cross  Section  of  a  Printed  Strip  Dipole  Excited 
by  a  Transmission  Line  Embedded  In  the  Dielectric 
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design  procedure  which  yields  the  microstrip  dipole  length, 
overlap  and  offset  so  that  a  desired  input  match  can  be 
achieved  for  a  given  substrate.  In  all  the  results  which 
are  going  to  be  presented  below,  the  operating  frequency  is 
10  GHz. 

A.  Self- Impedance  Evaluation 

For  the  evaluation  of  the  microstrip  dipole  self¬ 
impedance,  transmission  line  theory  is  applied  and,  because 
of  this,  unimodal  behavior  of  the  field  far  from  the  strip/ 
dipole  coupling  region  is  essential.  In  order  to  satisfy 
this  requirement,  the  transmission  line  is  kept  very  close 
to  the  ground  plane  giving  a  ration  W2/(h-bg)  >  2.0. 

For  the  particular  geometry  considered  here,  a  uni¬ 
modal  field  is  excited  under  the  transmission  line  and,  as 
a  result,  the  amplitude  of  the  current  distribution  beyond 
an  appropriate  reference  plane  looks  like  a  standing  wave 
which  is  due  to  TEM  waves  traveling  in  opposite  directions. 
For  this  reason,  the  microstrip  line  is  approximsted  by  an 
ideal  transmission  line  of  characteristic  impedance  which 
is  terminated  in  an  unknown  impedance  (see  Figure  5.3). 

This  Quasi -TEM  mode  has  a  wavenumber  B  and  a  standing  wave 
ratio  SWR  equal  to  the  average  values  evaluated  from  the 
original  current. 

If  the  origin  of  x  coordinates  is  taken  at  the  posi¬ 
tion  of  Z_,  then  the  voltage  and  current  waves  on  the  ideal 
transmission  line,  with  respect  to  this  plane  of  reference. 


«re  given  by; 


Figur*  9.3 

Current  Amplltud*  on  ttio  Strip  Dipolo  and 
Mierootrip  Tronomlooion  Uno 


-j6x  ^  g  jBx 


V(x)  ■  Ae 


I(x)  « 


(5.1) 


(5.2) 


If  one  considers  the  positions  and  of  two  con¬ 

secutive  maxima  and  minima,  equation  (5.1)  gives 


"max  .  "  - 


(5.J) 


Since  the  absolute  difference  between  x^^„  and  x_._  is 

max  inin 

equal  to  one-fourth  of  the  wavelength  of  the  guided  wave 
(x«e.v'X„._  •  Xg/4),  then  equation  (5.3)  can  be  written  as 

iii3x  iimi 


B  ^^®^max  ^.TT 

^  '  J  ^  e  ^ 

I  B  J2®^min 


The  reflection  coefficient  r  is  defined  as 


r(x)  -  I  e^^®*  -  r(0)e^^^* 


If  one  considers  that 

V. 


(5.4) 


(5.5) 


(5.6) 


then  equation  (5.3)  gives  the  following  expression  for  the 
reflection  coefficient. 


49 


1  ♦  r(x^) 


(5.8) 


If  now  this  plane  of  reference  is  considered  at  the  posi¬ 
tion  of  the  first  current  maximum  d  „  from  the  end  of  the 

max 

open  circuited  transmission  line,  then  the  self  impedance 
measured  with  respect  to  this  plane  is  given  by 


Z  fd  )  1  +  rfd  )  1  +  r[0)e 

s'-  max-'  _  max^  _  ^ 


1  -  r(0)e 


j2Bd. 


(5.9) 


Since  for  the  evaluation  of  r(d  ),  the  method  of  moments 

Hio  A 

together  with  transmission  line  theory  has  been  applied, 
equation  (5.9)  indicates  that  the  evaluated  self  impedance 
depends  on  the  characteristics  of  the  substrate  (e^,h)  as 
well  as  on  the  embedding  distance  of  the  transmission  line 
the  overlap,  the  offset  and  the  length  of  the  dipole. 

The  fundamental  design  procedure  is  now  revealed. 

One  wants  to  choose,  for  a  given  substrate,  the  right  posi 
tion  of  the  dipole  so  that  optimum  resonance  (in  other 
words,  perfect  match)  is  obtained.  The  condition  for 
optimum  resonance  is  characterized  by  the  relation 
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1 


(5.10) 


“■27- 

which  combined  with  (5.9)  gives 

r(d,„)  -  0  C5.li) 


or 

r(x)  -  0  (5.12) 

This  means  that  one  wishes  to  find  the  geometry  which 
perfectly  matches  the  dipole  to  the  transmission  line. 


B.  Dipole  Length  and  Overlap  Variation 
The  self  impedance  is  evaluated  using  equation  (5.9) 
as  a  function  of  the  length  of  the  dipole  and  the  over¬ 
lap  K, 


ovp 


ovp 


Si 
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(5.15) 


where  Sj  is  the  part  of  the  surface  S^j  of  the  dipole  which 
is  over  the  transmission  line  (Fig.  5.1).  The  real  and 
the  imaginary  parts  of  the  normalized  self  impedance  are 
plotted  (Fig.  5.4)  for  a  substrate  thickness  h  *  0.077", 
dielectric  constant  •  2.53  and  embedding  distance 
■  0.0485".  The  width  of  the  strips  is  w  ■  0.060"  and 
the  thickness  t  ■  0.00025".  The  length  of  the  dipole  is 
varying  between  0.373"  and  0.349"  and  the  overlap  takes 
values  between  38t  and  86t.  Figure  5.4  implies  that  as  the 
overlap  becomes  larger,  the  dipole  length  curves  approach 
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the  origin  of  the  axis,  which  means  that  as  the  dipole 

moves  to  the  left  from  the  end  of  the  transmission  line, 

the  coupling  weakens  and  the  reflection  coefficient  on  the 

microstrip  line  approaches  its  value  at  open  circuit. 

However,  there  exists  a  particular  overlap  for  which  the 

Z, 

curve  of  the  dipole  lengths  goes  through  the  point  ■  1; 
in  other  words,  for  the  geometry  considered  above,  there  is 
one  value  for  the  overlap  and  a  specific  dipole  length  so 
that  the  dipole  is  perfectly  matched  to  the  transmission 
line. 

C.  Dipole  Length  and  Offset  Variation 

It  is  of  interest  to  investigate  now  the  dependence 
of  the  normalized  self  impedance  on  the  offset  and  the 
dipole  length.  Again,  the  real  and  imaginary  parts  of  the 
normalized  self  impedance  are  plotted  (Figure  5.5)  for  a 
substrate  thickness  h  ■  0.077”,  dielectric  constant  = 
2.35  and  embedding  distance  b^  *  0.0485”.  The  width  of  the 
strips  is  w  ■  0.060”,  the  thickness  t  ■  0.00025”  and  the 
overlap  ■  501.  The  length  of  the  dipole  varies 

between  0.373”  and  0.349”  and  the  offset  takes  the  values 
0.00”,  0.02”,  0.04”,  and  0.08”.  From  Figure  5.6  one  can 
see  that  as  the  offset  becomes  larger,  the  coupling  weakens 
and  the  dipole  length  curves  are  shifted  to  smaller  values 
for  the  self  impedance.  Also,  the  resonant  length  of  the 
dipole,  or  in  other  words  the  dipole  length  for  which  the 
self  impedance  becomes  real,  is  a  monotonically  descending 
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0.041  \o 


Fuctlon  of  the  Offset  and  L 


of  the  offset,  which  means  that  as  the  offset  changes  to 
larger  values,  the  dipole  self  impedance  becomes  more 
capacitive.  Again,  as  happened  in  the  case  of  the  overlap 
variation,  one  can  see  that  there  exists  a  specific  offset 
which  can  give  optimum  resonance  if  the  dipole  length  is 
changed  appropriately. 

D.  Dipole  Length  and  Substrate  Thickness  Variation 

As  mentioned  before,  the  normalized  self  impedance 
is  also  a  function  of  the  substrate  characteristics  and 
the  position  of  the  transmission  line.  In  order  to  examine 
this  dependence,  the  width,  thickness,  overlap  and  the 
distance  of  the  transmission  line  from  the  ground  plane 
are  kept  constant  at  the  values  given  above,  the  dielectric 
constant  becomes  equal  to  2.53,  the  dipole  length  is  varied 
again  between  0.373"  and  0,349",  and  the  substrate  thick¬ 
ness  takes  the  values  0.0645",  0.0765",  0.0825",  and 
0.0945".  Figure  5.6  indicates  that  as  the  substrate  thick¬ 
ness  takes  larger  values,  the  coupling  between  the  dipole 
and  the  transmission  line  becomes  smaller,  the  dipole 
resonant  length  decays  monotonically  and  the  self  impedance 
for  a  given  dipole  length  becomes  more  capacitive.  Also, 
as  happened  in  the  two  previous  cases  of  the  overlap  and 
offset  variation,  there  is  a  unique  value  for  the  embedding 
distance  which  combines  with  the  appropriate  dipole  length 
to  give  a  perfect  match. 

It  is  also  interest'ng  to  see  how  the  current  dis- 
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tribution  changes  as  a  function  of  the  substrate  thickness 
h  around  optimum  resonance.  Figure  5.7  displays  the  cur¬ 
rent  distribution  normalized  to  the  incident  current  on  the 
transmission  line  for  h  ■  0.099"  (undercoupled),  h»0.093” 
(perfectly  matched)  and  h«0.082"  (overcoupled).  The  sub¬ 
strate  is  duroid  (e^  ■  2.53),  the  distance  of  the  trans¬ 
mission  line  to  the  ground  plane  is  0.0285"  and  the  length 
of  the  dipole  *  0.353",  As  shown  in  the  figure,  the 
current  on  the  dipole  takes  its  maximum  value  when  it  is 
perfectly  matched. 

E.  Comparison  with  Experimental  Results 

The  theoretical  analysis  of  a  printed  strip  dipole 
electromagnetically  coupled  to  an  embedded  microstrip  line 
is  tested  by  comparing  theoretical  results  to  experimental 
ones.  Stern  and  Elliott  (43]  measured  experimentally  the 
self  impedance  of  strip  dipoles  with  rounded  corners 
(Figure  5.8)  printed  on  duroid  boards  (c^  ■  2.35)  of  sub¬ 
strate  thickness  h*0.077",  excited  by  a  microstrip  trans¬ 
mission  line  in  the  dielectric  at  a  distance  from  the 
ground  plane  equal  to  0.0285".  The  width  of  the  strips  was 
w  ■  0,060”  and  the  thickness  t  ■  0.00025".  The  self 
impedance  was  measured  for  different  dipole  lengths  and 
its  normalized  values  are  shown  on  a  Smith  chart  (triangles) 
in  Figure  5.9.  The  solid  line  corresponds  to  the  theoreti¬ 
cal  results.  From  Figure  5.8  one  can  see  that  the  experi- 
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Current  Amplitude  on  the  Strip  Dipole  and  Trenemiselon  Line 


Figure  5.8 


Printed  Strip  Dipoie  with  Round  Corners 
EM  Coupied  to  a  Microstrip  T.L. 


mental  resonant  length  is  about  0.390"  while  the  theoreti¬ 
cal  one  is  0.379".  Therefore,  there  is  a  difference  of 
2.751  with  approximately  2%  resulting  from  the  different 
shapes  of  the  dipoles.  Those  which  were  studied  analytically 
h  '  a  rectangular  shape  whereas  those  measured  experi¬ 
mentally  had  round  corners  [56]. 

The  difference  between  the  theoretical  and  experi¬ 
mental  values  of  the  self  impedance  has  two  attributions: 
i)  Different  shapes  of  the  dipoles 
ii)  The  fact  that  in  the  theoretical  evaluation 
of  the  current,  the  hybrid  nature  of  the 
modes  propagating  in  the  microstrip  was 
taken  into  account  while  for  the  experi¬ 
ments  only  an  equivalent  TEM  mode  was 
measured. 

5-2.  Strip  Dipoles  Excited  by  a  Gap  Generator 

This  section  of  Chapter  5  presents  design  procedures 
for  microstrip  dipoles  printed  on  or  embedded  in  the 
dielectric  substrate.  The  dipoles  are  center-fed  by  an  in 
phase  unit  voltage  delta  gap  generator.  All  the  dimensions 
presented  are  normalized  with  respect  to  the  free  space 
wavelength  Due  to  an  assumed  time  dependency  of  e^“^, 

inductive  reactance  is  positive  in  all  plots.  The  dipole 
is  considered  either  alone  or  in  the  presence  of  parasitic 
dipoles  printed  or  embedded  in  the  dielectric.  The 
material  given  here  relates  the  antenna  geometry  (dipole 
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length,  substrate  thickness,  dipole-ground  plane  distance, 
position  and  length  of  the  parasitics)  to  antenna  charac¬ 
teristics  (resonant  length  and  resonant  resistance).  The 
presentation  of  the  numerical  results  is  completed  in  three 
steps:  At  first,  a  dipole  embedded  in  the  substrate  is 
considered  and  its  characteristics  are  discussed  in  terms 
of  the  embedding  distance.  After  that,  this  dipole  is  con¬ 
sidered  in  the  presence  of  one  printed  parasitic  dipole  and 
the  change  in  its  performance  is  studied.  Finally,  the 
same  dipole  is  considered  in  the  presence  of  two  parasitics, 
one  printed  on  the  interface,  the  other  embedded  in  the 
dielectric.  Its  characteristics  are  studied  in  terms  of 
the  relative  positions  of  the  parasitics  as  well  as  their 
overlap. 

A.  One  Dipole  Printed  or  Embedded  in  the  Substrate 

One  of  the  most  important  characteristics  of  a  dipole 
is  its  input  impedance  (Fig.  5. ID).  Figures  5.11  and  5.12 
show  the  real  and  the  imaginary  parts  of  the  input  imped¬ 
ance  when  the  strip  dipole  is  printed  on  a  duroid  board 
(Cy  *  2.45)  with  substrate  thickness  h  equal  to  0.2  and 
for  different  strip  widths  (w  *  0.0002X^,  O.OOIX^,  O.OIX^). 
Figure  5.11  shows  that  the  input  resistance  around  reson¬ 
ance  is  insensitive  to  width  variation  while  its  value  far 
from  resonance  becomes  lower  as  the  width  increases.  The 


effect  on  the  resonant  length  is  also  very  small  as  Figure 
5.12  reveals.  However,  the  input  reactance  is  very  sensi- 


(r«2.46,ht:0.2A.  and  t=  0.0001a 


Imaginary  Part  of  tlia  Input  Impodanoo  of  a  Printed  Strip 
Dipole  with  •r-2.46,h-0.2A.  end  t  •  0.0001a. 


tive  to  changes  in  the  width  talcing  much  lower  values  as 
the  width  becomes  larger.  From  the  figures  one  can  con¬ 
clude  that  the  behavior  of  the  input  impedance  of  a  printed 
strip  dipole  as  the  strip  width  is  varied  is  similar  to 
that  of  a  free-space  wire  dipole  as  its  radius  is  varied. 
Figure  5.13  shows  the  current  distribution  in  amplitude 
and  phase  on  a  strip  dipole  printed  on  duroid  (e^  ■  2.45) 
with  substrate  thickness  h  *  0.2X^,  w  *  O.OIX^, 
t  ■  O.OOOIX^  and  for  L  equal  to  0.62X  ,  0.347X  ,  and 

0  O  0 

0.248X  . 

0 

Consider  now  the  dipole  shown  in  Figure  5.10  for 
duroid  of  dielectric  constant  ■  2.53,  substrate  thick¬ 
ness  h  ■  0.065X  ,  thickness  of  strip  t  ■  O.OOOIX  and 
o  o 

width  w  ■  0.05X^.  Figure  5.14  shows  how  the  resonant 
length  Ly  and  the  resonant  resistance  vary  as  a  function 
of  the  embedding  distance  b^.  It  is  interesting  to  note 
that  as  the  strip  dipole  enters  the  dielectric  substrate 
and  moves  closer  to  the  ground  plane,  the  resonant  length 
decreases  to  a  minimum  value  when  the  dipole  is  at  an 
embedding  distance  equal  to  the  half  of  the  substrate 
thickness,  and  after  that  increases  to  a  maximum  value 
when  the  dipole  approaches  the  ground  plane.  However, 
the  behavior  of  the  resonant  resistance  is  different.  It 
takes  its  maximum  value  when  the  dipole  is  printed  on  the 
interface  and  decreases  to  zero  as  the  dipole  enters  the 
interface  and  moves  close  to  the  ground  plane.  The  reson- 
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Figure  6.14 

R«aonant  Langth  and  Raaonant  Raaiatanca 
of  a  Strip  Dipola  aa  Functiona  of  tha  Embaddlng 
Diatanea.  <r  ■2.6S,li  s  O.OSSa.  .t«  0.0001a. 
and  w>  0.06a. 


ant  resistance  when  the  dipole  is  on  the  interface  takes 
such  a  small  value  because  the  substrate  is  thin  and  the 
strip  dipole  is  wide. 

B.  One  Embedded  Strip  Dipole  Excited  by  a  Gap 

Generator  in  the  Presence  of  a  Printed  t>arasitic 
bipole 

The  dipole  of  Figure  5.15  is  now  considered.  In 
this  case,  the  relative  dielectric  constant  for  the  sub¬ 
strate  is  2.35,  the  substrate  thickness  h  is  equal  to 
0.065Xp,  the  strip  thicknesses  are  ■  t2  •  0.0001,  the 
strip  widths  are  w^^  *  w^  ■  0.05X^  and  the  overlap  is  23%. 

As  shown  in  Figure  5.16,  the  resonant  length  of  the 
excited  dipole  as  a  function  of  the  embedding  distance 
does  not  go  through  a  minimum  as  happened  in  the  case  of 
the  single  dipole  but  decreases  monotonically  as  the 
exciter  approaches  the  ground  plane.  The  behavior  of  the 
resonant  resistance  is  not  that  much  different.  For  values 
of  b^  smaller  than  half  of  the  substrate  thickness,  it 
oscillates  around  a  value  of  15  ohms  and  for  b^  larger 
than  h/2  it  decreases  monotonically  to  zero  as  the  exciter 
goes  very  close  to  the  ground  plane.  It  is  interesting 
to  compare  the  variation  of  the  resonant  length  and 
resonant  resistance  as  functions  of  the  embedding  distance 
for  the  cases  of  the  single  dipole  and  the  dipole  in  the 
presence  of  a  parasite  (Figures  5.14,  5.16).  One  can 
observe  that  the  parasitic  not  only  changes  the  behavior 
of  the  resonant  length  and  resonant  resistance,  but  also 
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Figure  5.15 


Embedded  Strip  Dipole  Excited  by  a  Voltage 
Gap  Generator  In  the  Presence  of  a  Parasitic 
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Figure  5.16 

Resonant  Length  and  Resonant  Resistance 
of  the  Excited  Dipole  as  Functions  of  its 
Embedding  Distance  Cra2.53  ,  h  «0.065x. 
0.06 X.  1^23%  and  t-t-  0.0001  x. 


makes  their  values  larger.  This  means  that  the  band¬ 
width  and  the  efficiency  of  the  antenna  are  increased.  It 
is  also  interesting  to  see  how  the  resonant  length  and 
resonant  resistance  change  as  functions  of  the  substrate 
thickness  when  the  distance  of  the  exciter  from  the  ground 
plane  is  held  constant.  In  Figures  5.17  and  5.18,  the 
resonant  length  and  resonant  resistance  are  shown  as 
functions  of  h  when  Cj.  ■  2.35,  w^  ■  W2  ■  0.05X^,  *  t2  “ 

O.OOOlX^and  h-bg  ■  O.IX^,  0.2X^  respectively.  In  the 
first  of  these  two  figures,  the  substrate  thickness  takes 

values  between  0.15X  to  0.25X  while  in  the  second  one  h 

0  o 

varies  between  0.065X^  and  0.075X  .  Except  for  the  differ- 

0  0 

ence  in  values  the  behavior  of  L_,  R_  is  the  same  in  both 
cases.  From  Figures  5.17  and  5.18  one  can  see  that  as  the 
parasitic  dipole  moves  to  greater  distances  from  the  exciter 
the  coupling  weakens  and  the  resonant  length  and  resonant 
resistance  of  the  embedded  dipole  asymptotically  tend, as  one 
would  expect  to  the  values  they  have  in  the  case  of  the 
single  excited  dipole  (Figure  5.14). 


C.  M  Embedded  Strip  Dipole  Excited  by  a  Gap 
Generator  in  the  Presence  of  Two  Parasitic 
Dipoles 

The  geometry  for  this  antenna  is  shown  in  Figure  5.19 
For  this  arrangement,  with  ej,  ■  2.35,  h«  0.065X^,  bg  • 
0.04U^,  tj  -  .  O.OOOU^kW  -  -  sot,  Wj  -  Wj  - 

O.OSX^f  the  real  and  imaginary  parts  of  the  input  impedance 
for  two  different  values  of  8  (distance  of  the  embedded 
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Figure  6.18 

Resonant  Length  end  Reeonent  Reeletence  of 
the  Excited  Dipole  ee  Funetlone  of  the  Substrate 
Thlekncas  c,.2.3S.  h-b,-  0.02Sx.  ,  w  >  O.OSx. 
•c;^23%  and  Vt..  0.000 n. 
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Figure  5.19 

Embedded  Strip  Dipole  Excited  by  a  Voltage  Gap 
Generator  In  the  Presence  of  Two  Parasitica.  One 
Embedded  in  the  Dielectric  and  the  Other  Printed 

on  the  Interface 
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parasitic  from  the  interface)  are  given  in  Figures  5.20  and 
5.21.  From  these  two  figures,  one  can  see  that,  as  the 
parasitic  approaches  the  exciter,  the  resonance  occurs  at 
longer  lengths  of  the  excited  dipole,  while  the  bandwidth 
decreases  and  the  resonant  resistance  increases  to  values 
above  100  ohms.  The  resonant  length  and  resonant  resistance 
as  functions  of  the  embedding  distance  6  are  shown  in 
Figure  5.22  for  two  different  values  of  the  overlap 
(501,  80t). 
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Figur*  S.20 

Input  Imptdmnc*  of  tho  Excitor  no  o  Funotion  of  Ito 
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CHAPTER  6 


RADIATION  PATTERN  OF  PRINTED  STRIP  DIPOLES 

6*1.  Formulation  of  the  Equations  for  the  Far~Zone  Field 
Radiation  pattern  or  far  zone  fields  can  be  obtained 
by  using  a  rigorous  numerical  technique  as  described  in 
Chapter  2.  The  numerical  technique,  although  valid  for  all 
distances,  is  quite  expensive  for  the  far  zone  computations 
In  the  calculation  of  the  radiation  pattern,  only  the  far 
zone  fields  above  the  air-dielectric  interface  are  required 
Hence,  the  conventional  stationary  phase  method  I57]-[58] 
can  be  used.  The  details  of  the  method  are  covered  in 


Appendix  D.  When  a  printed  strip  dipole  is  excited  either 
by  a  gap  generator  or  by  a  microstrip  line  embedded  in 
the  substrate  with  h  -  b_  <<  X  (see  Figure  6.1),  the  far- 
zone  field  is  due  totally  to  the  radiation  from  the  dipole. 
Under  this  assumption,  the  far-zone  electric  field  is  given 


t(R,e,«) 


-jk.P 


(6.1) 


In  equation  (6.1)  R_  is  a  scalar  (called  the  antenna  shape 

factor)  and  is  given  by 

cos(k,l  sinecos4i)  -  cos(kt  ) 

Ra  -j260  - LJS - —11^ - 2^  . 

k  sin(ki„)  1  -  sin^ecos^^ 


(kj  ^  sin^sinej-  ~  S(wJ.e,^) 
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Figure  6.1 

Printed  Strip  Dipole  Excited  by  a  Microetrip 
Tranemieeipn  Line  Embedded  in  the  Dielectric 


jkjajSinecos4i 

*  e  e 


n  jk,  (n£„)sinecos4> 


(6.2) 


where  is  the  zeroth  order  Bessel  function  and  S(w^,e,4)) 
o  e 


is  the  integral 

S(Wg,e,<t»)  -  j 


rcos'^Cw/wM  Wg 

cosfkj^  -j—  sin^sinScosaJda 

'0 

(6.3) 


Also,  is  a  vector  called  the  substrate  factor  and  has  the 


form 


*5  ■  ^''se  *  K* 


(6.4) 


where 


Rgg  »  -coS(|)  ♦(e^,h,0)  tcosO  ♦  (Cy-l)sinetaneA  (c^,h,e) ) 

(6.5) 


Rg^  -  sin(>  ♦(e^.h.e) 


(6.6) 


and 


♦(Cy.h.e)  -  ^  . 


cose 


I  COS0- j /ej.“Sin^e  cot  (kj/ej.-sin^e  h) 


(6.7) 


A(ej.,h,0) 


cose 


:^cos0+  j»4y*sin^e  tan(kj^/cy-sin^e  h) 


(6.8) 
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6-2.  Effect  of  Substrate  Thickness  and  Permittivity  on 

the  Radiation  Pattern 

An  investigation  of  the  expressions  given  for  the  far- 
zone  electric  field  indicates  that  the  effect  of  the  sub¬ 
strate  properties  on  the  radiation  pattern  is  controlled  by 
the  factors  ♦(Ey,h,e)  and  A(ej.,h,e).  Furthermore,  it  is 
verified  that  ♦(e^,h,e)  is  a  result  of  the  substrate  guided 
TE  modes,  while  A(ey,h,e)  is  due  to  the  TM  modes.  A 
thorough  analysis  of  *  indicates  that  it  determines  the 
position  of  the  nulls  and  the  principal  as  well  as  secondary 
maxima  of  the  pattern  for  6  <  y.  If  one  considers  the 
radiation  pattern  of  a  strip  dipole  at  resonance,  then  the 
number  of  the  lobes  and  the  position  of  the  nulls  are 
totally  controlled  by  the  substrate. 

A.  Number  of  Lobes 

The  number  of  lobes  in  the  radiation  pattern  can  be 
determined  from  the  equation 


where  the  brackets  indicate  the  integer  value  of,  N  is  an 
integer,  and  a  is  an  arbitrary  real  positive  number  less 
than  one.  From  equation  (6.7),  the  following  cases  can  be 
identified: 

1)  If  N  •  0,  there  exists  a  single  lobe  with 
maximum  at  6  *  0. 

2)  If  N>0  and  a>0,  there  exist  2N  +  1  lobes 
with  one  of  the  maxima  at  6  ■  0. 
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3)  If  N>0  and  a'»0+,  there  exist  2N  lobes 


with  a  null  at  6  *  0. 

B.  Positions  of  Nulls  (N  >  0) 

From  equation  (6.9)  one  can  derive  the  relation 

If  it  is  assumed  that 

*  ” 

then  the  position  of  the  n^^  pair  of  nulls  is  given  by 

®n  1,2  *  y^r 

where  <  *  m+n  and  n  ■  1,2,...,N. 

The  dependence  of  the  f  and  fi-plane  normalized  radia¬ 
tion  power  patterns  on  and  h  can  now  be  investigated. 

For  a  duroid  substrate  with  Cj.  ■  2.35  and  for  a  substrate 
thickness  of  h  •  0.21^,  equation  (6.9)  is  satisfied  for  N®0 
and  a  ■  0.613  and  therefore  the  radiation  pattern  consists 
of  a  single  lobe  as  shown  in  Figure  6.2a.  In  light  of 
(6.7),  it  can  be  verified  that  for  h  ■  0.2Xj^,  0.9751^  and 
l.OSX^  the  J  and  plane  normalized  power  patterns  will  have 
one,  two,  and  three  lobes  respectively,  as  shown  in  Figures 
6.2b,  6.3a,  6.3b. 

If  the  substrate  thickness  h  is  fixed,  e.g.,  at 
h  ■  0.1016Xp,  then  the  f  and  fi-plane  normalized  power  pat¬ 
terns  are  shown  in  Figures  6.4a,  6.4b  and  6.5a  for  *  2, 


h 

TT 


\ 


*  N 


(6.10) 


(a)  h-hopt  h=0.2x 


10,  35  respectively.  With  increasing  e^.,  we  observe  that 
the  PCD  directivity  is  decreased  because  more  energy  is 
radiated  close  to  6  •  ^  direction  along  the  length  of  the 
antenna  as  the  number  of  nodes  guided  in  the  substrate 
increases.  It  is  further  observed  in  Figure  6.5b  that 
when  Cy  ■  25  and  h  ■  0.10161^  ,  there  exist  three  lobes  and 
according  to  (6.11)  the  nulls  are  0  »  ±62.111*.  This  last 
case  has  uncovered  an  interesting  phenomenon  which  is  the 
radiation  at  angles  6  very  close  to  y  .  This  will  be 
described  in  the  next  section. 

6-3.  Radiation  along  the  Horizon 

As  mentioned  above,  an  interesting  phenomenon  is  the 
existence  of  nonzero  radiation  close  to  the  horizon  in 
either  the  or  the  I  plane  when  the  electrical  thickness  of 
the  substrate  satisfies  specific  criteria.  In  particular, 
^-plane  radiation  along  the  horizon  is  seen  to  occur  when 
a  TE  surface  mode  turns  on  in  the  substrate  and  ^-plane 
radiation  close  to  the  horizon  is  seen  to  occur  when  a  TM 
mode  is  excited  and  the  substrate  can  support  more  than  one 
■^be.  This  phenomenon  is  explained  analytically  by  the 
coincidence  of  a  pole  and  a  branch  point  in  a  Sommerfeld- 
type  integration. 

'  In  order  to  highlight  the  criteria  under  which  radia¬ 
tion  along  the  horizon  takes  place,  the  case  of  the  dipole 
of  Figure  6.1  will  be  considered  in  detail.  If  one  studies 
equations  (6.1)  -  (6.8),  it  can  be  seen  that  the  ^and 


plane  radiated  fields  are  given  by 

-jkjR 


Eg(t«plane)  ■ 
*  jccse  +  (Ey- 


RgC^'O).  h,e)  • 


Sin^e 

cole- 


r.h,e)j 


E^(]?-plane)  ■  - 


-jklR 


“RaC<>  “  y)  ♦(ej.,h»6) 


From  equations  (6.12)  and  (6.13),  one  can  obtain 


•jkjR 


^a^^  “  7  »  ^  “  7^  * 


(6.12) 


(6.13) 


*^7  /Ej.-sin^e  cot (k2»4j.“Sin^e  h) 


(6.14) 


In  words,  the  radiated  field  in  the  H-plane  tends  to  zero 
at  the  horizon  (e-^Tr/2)  unless  the  denominator  also  tends  to 
zero  as  ©-►7r/2.  When  cot(kj/t^^  h)  is  equal  to  zero  then. 


2v/JrpT  -  CZm^  1)  j 


(6.15) 


^  «  iBjl  m  •  0,1,2,... 


(6.16) 


Equation  (6.16)  is  simply  the  condition  for  a  TE  surface 
wave  mode  turning  on.  Equations  (6.14)  and  (6.16)  imply 
that  there  is  nonzero  radiation  at  the  horizon  in  the  if* 
plane  pattern  when  a  TE  surface  wave  mode  turns  on.  For  the 
{•plane  pattern  one  finds  that 
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-jk  R 
e  *■ 
~k“ir" 


Ra(0- 


0,  e-J)  (Ej.-!) 


lim 


cose 

TVTT 


«  0 
(6.17) 


Therefore,  at  6  ■  -y  the  radiation  in  the  E-plane  is  always 
zero.  However,  it  has  been  observed  that  for  6  *  7  ■  b, 
where  b  is  a  very  small  angle,  the  radiation  is  nonzero  when 


— r  n  .  m 

0 


m  ®  1,2, .. .  (6.18) 

and  when  the  substrate  can  support  more  than  one  lobe,  i.e.. 


o 


h 

oj 


+  a 


(6.19) 


Equation  (6.18)  combined  with  (6.19)  gives 


♦^Cy- 1  1,2,...  (6.20) 

Equations  (6.16)  and  (6.20)  provide  a  set  of  curves  (see 

Figure  6.6)  which  one  can  use  in  order  to  choose  the  right 

substrate  so  that  the  dipole  radiates  close  to  the  horizon 

either  in  the  ?-plane  or  in  the  fi-plane.  Figure  6.7  shows 

the  H-plane  pattern  of  a  strip  dipole  at  resonance  printed 

on  a  substrate  with  c^.  ■  2.1  and  h  «  0.238X^  and  Figure  6.8 

shows  the  if-plane  pattern  of  the  same  dipole  printed  on  a 

substrate  board  with  e.  ■  2.286  and  h  ■  0.6614X^.  Also, 

r  o 

Figure  6.9  shows  the  E-plane  pattern  when  Cj.  *  4  (Quartz) 

and  h  ■  0.285X  . 

0 


91 


-90 


H 


X' 


Figure  6.7 

plane  Radiation  Pattern  of  a  Printed  Strip  Dipole 
with  €r.2.1.  h*  0.238  X*  and  w.t«0.000lA. 
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o 

o 


Figure  6.8 

H  -  plane  Radiation  Pattern  of  a  Printed  Strip  Dipole 
with  Cr>2.266,  ha0.661Xo  and  w«t»0.000lAo 
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Figure  6.9 

E  -plane  Radiation  Pattern  of  a  Printed  Strip  Dipole 


with  €r*4  ,  hm0.285Ae  and  w«tcO.OOOlA. 
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APPENDIX  A 


FORMULATION  OF  POCKLINGTON’ S  INTEGRAL  EQUATION 


2.35)  IN  THE  FORM  GIVEN  BY  (2.37 


In  Chapter  2  it  was  shown  that  the  integral  equation 
for  the  electric  field  is  given  by 


Ei(?) 


rw''/2  fK 

E  /  dy*  J  dx‘ 
v*l,2  J'V  /2  "^o 

r  2pi  .2pi  "I 

,.2,i  .  ^  *^vxx  .  ®  *^V2X| 

^iNxx  *  ^-7'  * 


J^(x*.y»,z'').  (A.l) 


If  one  considers  the  relationship 


3F^ 

vzx 


dx 


(A.  2) 


then  equation  (A.l)  can  be  written  as 

E^Cr)  ■  E  I  dy’  I  dx* 
*  vl,2  J-wV/2  Jo 


(A.  3) 


From  equation  (A. 2)  the  function  F^^  is  given  by 

■  -/->*  It  Kzi  ■  ■/"*  Kzi 

By  substituting  (2.28)  into  (A. 4)  one  can  obtain 

/iuy  \  f 

-  2(^j(e^-l)  fix  cos* 

>  If.  „  «  -u„i*«u5inh[u(h-i'’*i'’«^j)) 

TPT”  1 1  J-i  lAPje  — T—n — -  -v'T"' 

’i  l-'o 
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COSlu(h-Z  +2 


.  2 » Ey I 


r  r®  6., 

2^— Jdx  cos^J  Jj(Xp)e  ° 

|6ijUoCOsh(uh)  -  H 

fl(A.e^.h) 

[fi^jSinhCuh)  ♦  6^2®^*'^!'^ 3  2 


From  the  relations 


p  «  [(x-x’)2  ♦  (y-y’)^l’* 


X‘dX  (A. 5) 


(A,  6: 


cos4> 


it  is  determined  that 


-cos4» 


(A.  7] 


(A.S: 


dJjjCXp) 

a(Xp3 


1  dx’ 

X  HIP  dp 

1  1 

X  djp  cos4i 


-  -Jj^p) 

A  aubstitution  of  (A. 8)  and  (A. 9)  into  (A. 5)  yields 
FL(T/r')  •  -1)  f  J  (Xp)e 


(A.  9; 
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APPENDIX  B 


EVALUATION  OF  THE  SPACE  INTEGRALS  IN  THE  INTERVAL 

XetO.A] 


As  shown  in  Chapter  4,  the  elements  of  the  generalized 


impedance  matrix  are  given  by 


r  r  fA"''CX,e^.h;z^2'') 

Jo  v-1,2  f^C^,Ey,h) •f2tX,E^,h) 


mn  o 


where 


CB.l) 


{j  (Xp)}  -  iiZizii -  r  ^  dy'  J^(y’)  • 

i  °  Isin(kt^)]^  J-w'’/2  ^ 


*  dx  r  ^  dx*  sinlk(£  'X)l  sinlkCtj^-x* )  1  • 

Jck 


ISCp-Oj)  (P"P2)+5 (p-Pj)  +  SCp-p^)!  ♦ 


+  E2  *  dx  sinlk(tj^-x)]  I6(x’  +  tj^)  +  6(x»-£j^)  - 

-2cos(kS.^)6(x')]  .  I5(p-pj)  +  6(p-pj)lj  •  J^^CXp) 

(B.2) 

with  PjjP^rPjfP^  given  by  equations  (4.13)  (4.16)  and  p 

given  by 

p  -  {(x-x'+x^-x^)2  ♦  (y-y’)^}**  (B.3) 

The  Bessel  function  in  equation  (B.2)y  because  of  (B.3),  can 
be  written  as 
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i  / 

•^r\ 


TT  jX  (x-X’+Xjj^-Xj^)COSe 


COS  (X (y-y' )sine)de 
(B.4) 


If  one  substitutes  (B.4)  into  (B.2)  and  interchanges  the 
order  of  integration,  the  multiple  space  integrals 

V 

{J^(Xp)}  take  the  form 


mn 


L  {J-(xp)} 
mn 


«Cy-y‘)  .  1  r^e 

Jc\ 


(sinCkt^)]"^  ^  -D 


J  ^  dx  J  ^  dx’  sinlk(£jj-x)  ]sin[k(£jj-x' )  1  » 


jgjX(x+x')cose^  ^jX (x-x’)cose^  gjX(-x+x’)cose  ^ 

^  ^jX(-x-x*)cosej  ^  ^-jXi^cose  _ 

-  2cos(k£j^)j  •  dx  sinlk(£j^-x)] 

+  •  I  dy’J^(y’)cos(X(y-y’)sine) 

CB.5) 

Equation  (B.5),  with  the  use  of  the  equality 

2cosa  ■  e^®  +  e”^® 


can  give 


■;  .  J  Pde 

■n  ®  Isin(kl„)K  *4) 


(B.6) 


4  f\.  ^3X(x^-x„)cose 
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I 


v^/2 

v^/2 


dy'  Jy(y')  cos(X(y-y')sine) 


/: 

K  r 


dx  sinI>(£jj-x)]cos(Xxcose) 


dx*  sinIl»(£jj-x’)lcos(Xx'cose)  ♦ 


♦  E2lcos(Xfcj^cose)  -  cos(ktjj)l 


(B.7) 


In  equation  (B.7)  if  one  applies  straightforward  integra¬ 
tion,  interchanges  the  order  of  integration  and  assumes  that 


rw'’/2 


dyj:;(y)Joap,„)  - 


/2  >- 

'***'’  *mn  *  then 

equation  (B.7)  can  be  written  as 


(B.8) 


■  Jb|*2p(T;) 


t„\2<J  j2a 

— 2-_  J  (Xx  )  ♦ 
jyp?  o'  mn' 


♦  B. 


(it)  TT^  *^^0^^’Sn(nn)''^0'^W 


*m(n-l)^^ 


(B.9) 


where  A2p»  B^^  are  constant  coefficients. 
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APPENDIX  C 


FORMULATION  OF  THE  INTEGRAL  (4.20)  FOR  THE  TAIL 

CONTRIBUTION 

If  we  choose  A  in  such  a  way  that 

cothr(A^-k2)’*hl  *  1 

then  in  the  integral 


(C.l) 


mn'^ 


with 


JA  v-1,2  f 


L  Uo(xp)> 
mn 

(C.2) 

A^'’CX,c^.h;2\2'’)  -  tEj(]c2-le2)^E^,,j./'^^^(X,E^,h;2\2'’) 


7  ^  i 

+  [Ejk^-E,!*]  (X,c  ,h;z\z'’ 

V2 


(C.3) 


and 


r  (X.e,.h;z'.zV)  -  zfe)  . 

vxx  ^  VAirkf/ 

•  *  UpSinh(u2'’)  1 }  • 

r'  (A.c^.hjzi.z'’)  .  2(^) 

V2  ^  Wkt/ 

(fiijUoCOshCuh)  -  d * 


sinhlu(h-z''  ♦  1  * 

/j“Pn\  r  1. 

e  ”  (Ej-l)  • 


{fi^jSinhCuh)  ^ 


(C.5) 
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the  functions  fj(^,ej.,h),  f2(X,ey,h)  can  be  approximated  by 
fj(X,Ey,h)  -  sinh(uh)  CUjj+u)  CC.6) 


“  sinh(uh)  (ej.u^+u) 


(C.7) 


In  equations  (C.2)  (C.5),  the  z-coordinate  for  the  source 

currents  (z^)  and  observation  points  (z^)  is  given  by 


t  i-1 

O 

€ 

n 

-bg+t  i-2 

z'*  »  ( 

-b.  v=  2 

i"  «  \  ^  ‘  {  (C.8) 

1.2j  |.b3  V.  2| 

If  one  substitutes  (C.6)  and  (C.7)  into  (C.4),  then  the 
function /^^j^{X,ey,h;z^,z'’)  can  be  written  as 


vxx  V4iTkt/ 1  “ “o  , 


o  1  ♦  H-> 
% 


-Upt. -ubg 

*  (»il«v2  *  «i2«vl’5  -  * 

1 

.  ,  ,  X  -ut  .  ^ 

1 2  v2  u  u 

1+  — 
u 


Similarly,  the  function  (X,c  ,h;z^,z'’)  can  be  approxi- 

vz 

mated  by 

^  i  2  .. 


j 


.Ej.,h) 

(C.9) 


r  -  2{^)  . 

vz  ^  \4iTkr/ 
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■*  («il«v2 


I- "'’7-^ - ‘ 

O  Vl*  —  E  ♦  i!_  J 

'  “o  ’■  “o/ 

.  ,  ,  ,x  .-V.-“'>s/  1  1 

*  *i2«vlJS  *  *  1 - u- - u-l* 

\1  ♦  -2.  1  +6  — / 

\  u  ^  *'r  U  / 


♦  «12'v2S* 


(■  -V  ■  ■  -Xi) 


(C.IO) 


In  equations  (C.9)  and  (C.IO)  u^,u  are  functions  of  X  given 
by 

(C.ll) 

(C.12) 


“o  * 


u  «  (X^  -  k^)** 


Because  of  (C.IO)  and  (C.ll),  the  functions  1  » 

u 

1  ♦  ^  ^  »  1  ♦  E  ~  in  (C.9)  can  be  written  as 

“o  o  *  “ 


1  +  ^  •  4  {1  -  ©^(x)}*^ 

(C.13) 

o 

1  ♦  •  Y  {1  ♦  e4(X)}’^ 

(C.14) 

'r  *  S-  *  nVr  <1  - 

(C.15) 

0  r 

1  *  'r^-  rVr  <1  * 

(C.16) 

where  ,  , 

k^  -  k^ 

0<e2(X)<J  \  (C.17) 

A  * 


409 


“  i  i  r(s;*D  • 


K2  Kj 


k^-k^ 

1  ^1 

0  <  e^CX)  <  j  •  2  "2 

A  -k2 

2e  k^-kj 
1  i  e,(x)  <  r-rr;  • 


(C.18) 

(C.19) 

(C.20) 


Also,  in  equations  (C.9]  and  (C.IO),  the  function 
■  (Uot  ♦  ubg)  approximated  by 


-(u„t  ♦  ubp  J. 

e  *  e 


(C.21) 


where 

t*  -  bg  ♦  tIl+2e4(X)  1  (C.22) 

If  one  substitutes  (C.13)  -  (C.16)  and  (C.21)  into  (C.9)  and 
(C.IO),  equation  (C.2)  takes  the  form 


■  I 


/, 


7  7  1*^1^^^® 


-“o* 


4xkj/  v*l,; 

V  r"» 

^  {J  (Xp)}du  *  I  2  2  k 


-ut< 


• ^  {J«CXp))  du  ♦ 


mn 


£  {Jo(J^P)>  du 
mn 


,  *u(2b  -t) 

;(X)e  s  . 


(C.23) 
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where 


Dj(X)  -  jlEj(kJ->2)+E2lel  +  lE2lf^-E2Vl  *  fiijfivl’ 

(— - ^ - 

\l-e2(A)  Cj.-1  l-ejCA)/! 


(C.24) 


D2(X)  -  ItEjCkJ-k^)  ^  £^,,3 


^i2*v: 


1  ♦  e.(X) 


♦  tEiV"-E2VKfiii6^2*«i2^1^  ‘ 
■(l*e4(X)  ■  c;^  •  r*e^irj)| 


(C.25) 


DjCX)  -  -  lEjCk^-V^)  ♦  E2le] 


l^e^CX) 


♦  64,6 


i2‘'v2 


(l‘eju)  '  1^  •  l*ej(X))j 


(C.26) 


From  all  the  equations  above  one  can  conclude  that 


|Di(X)-Dj(A)j  <  |Dj(«)-Dj(A)|  -  OCA’2) 


(C.27) 


Therefore,  equation  (C.25)  can  be  written  as 


4  /lup  \  f  “U  t  V 

Mirk*/  v-1,2  I  ^  J(A^-kJ)**  mn  ®  ®  ^  « 
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*  ®2(*)  I'l  2  •>  ‘'“‘V  * 

^  JCA^-k;)**  inn 


-u(2b_-t) 

(A‘-k2)'  ron 


where 


•  “5^*^  /(A^-kl)’’  ' 

p*  •  p(l  ♦  ej(A)) 


(C.28) 


(C.29) 


and 


p*  -  p^l  ♦  ej(A)J 


(C.30) 
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APPENDIX  D 


■FAR-ZONE  ELECTRIC  FIELD  OF  A  PRINTED  STRIP  DIPOLE 

The  far-field  due  to  a  printed  strip  dipole  (see 
Figure  6.1)  is  given  by 


and 


where 


and 


Eq 'V' kjlcosecos^n^  -  sinen^l 


E^  kj[-sin((illj^l 


(D.l) 


(D.2) 


jitiU  r  L+a, 

n,  - - %  /  ^dx*  /  ^  dy*  J(x’.y’)  • 

*  4iik^  -'-w/24bj 

(D.3) 


n. 


/L+a^  rw/2+b, 

■‘dx’/  dy' 


J(x’,y’)  /  h{^^(Xp)  e 


.  sinh(uh) 


cosh(uh) 


A^dX 


(D.4) 


In  equation  (D.3)  and  (D.4),  J(x’,y')  is  the  current 
density  on  the  strip  dipole 
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with  given  by  equation  (3.21)  and  fj(X,ey,h), 

f2(3^,ey,h)  given  by  equations  (2.31) »  (2.32).  If  the  sub¬ 
stitutions  X  -  k^sina,  z  ■  rcose,  p  •  rsine  are  introduced 
(where  r,e,^  describe  the  far-field  in  spherical  coordin¬ 
ates)  together  with  the  large  argument  asymptotic  expansions 
for  H^^^(Xp),  Hj^^(Xp),  Equations  (D.3)  and  (D.4)  can  be 


written  as 


iuv  rrr 

n  «  .  /-^  I  I  dy*J(x*,y’)  • 

^  yxrkj^  •'-w/2+bj 


JT^*  fCe^.h.a) 


- jrkjcos(e-a) 


(D.6) 


/•  L'^'a.  ^w/2^b. 

.  sid-ocos*/  dx'/  V' 

Tx  ywrkj  r  •'-w/2+bj 

/♦OB  -  jrk.cos  (e-a) 

g(Cj.,h,a)  e  da 

m 


(D.7) 


where 


f(ey.h,a) 


sinacosa 


sin(kj/ej.-sin^a  h) 


^sindsina 


and 


gC^y.h.a)  -  -k. 


.  2 

sin  a  cosa 


'sinesina 


sin(kj»4j.“sin^a  h) 
fj(a,Ey,h) 


cos (kj/Ey-sin^a  h) 

f2(a,ej.,h) 


CD.  9) 


The  integrals  with  respect  to  a  in  equations  (D.6)  and  (D.7) 
have  a  stationary  point  at  a  ■  6.  Therefore,  one  can  apply 
the  stationary  phase  technique  for  the  evaluation  of  these 
integrals  and  are  given  by 

_  rL*a.  -w/2+b, 

rrr"  I  ^  f  ^ 

dx’J  dy'J(x',y’)  • 

X  4x  prkj  aj  -'-w/2+bj 


.  f(e 


j,h,e)  j 


♦®  “ jrkj^cos (e-a) 


(D.IO) 


n  -  -  i/ili~  (l-ej.)cos^  J  dx*  f  dy’ 

z  ^47r  ^vrkj  r  -^-w/Z+bj 


•w/2+b. 


-jrk,cos  (6-a) 

•J(x',y')  g(Ej.,h,e)  e  da  (D.ll) 


1  .  1 


If  the  substitutions  r  ■  R-x'sinScosiJ'-y’sin^isine,  then 

the  stationary  phase  method  yields  (D.12) 

"x  -  -  TiT  hjT-  •  I 


Bj  •  Sj_u —  (Cj.-ljcosnane  «(c,tliie). 

1  .  ^  A  N  » 


A  C^r  » ® ) *  I 


(D.14) 
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where 


cose 


♦(Ej-.^tS)  *  ^  ‘  I  _  _ 

^  ^cose-j/cy-sin^e  cot(kj^/ey-sin^e  h) 


A(Ej,h.e) 


(P.15) 


cose 


CyCose+j v4j.-sin^e  tanCkj/Cj.  -sin^  h) 


(D.16) 


and 


jk,x*  sinecos4»f'^^^*^l  jk,y*sin(fisine 
dx'  e  ^  I  dy‘  e  ^ 

•'a  •'“W/2+b, 


N 


'W  “ 


(0.17) 


After  straightforward  integration  the  integral  I  takes  the 
form 


I  -  2 


cos  (kjtj^sinecos4>)-cos 


Ifsin(ktjj) 


2  2 
1-  —4  sin  ecos  ^ 

L 


►  • 


•|jo(l«i  ^  sin^sine)-  I  S(wJ.e.^) 

jk^b^sin^sine  jk^^a^sinecosi^ 
c  © 


n  jk, (n£„)sinecos^ 

£ 

n- 1  " 


where 


1  fcos‘^(w/w^)  r  wj 

SCwJ.e.^)  -  J  *  cos[kj  ^ 


(D.18) 


sin^sindcos 


ol  do 
J  (D. 


19) 
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(D.20) 


From  all  the  above,  one  can  conclude  that 

-jk^R 

where 

R^  -  j  60  I  (D.21) 

and 

K  -  «s6  * 

with 

R^e  «  -cos^  ♦(£y,h,6)rcos6  ♦  (E^-l)sinetane  • 

•  A(Cj.,h,6)  1  (D.23) 

and 

Rg^  «  sin^  •(Cy,h,e)  (D.24) 
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